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Abstract. We study the cosmology of bimetric theory with a composite matter coupling. 
We hnd two possible branches of background evolution. We investigate the question of stability 
of cosmological perturbations. For the tensor and vector perturbations, we derive conditions 
on the absence of ghost and gradient instabilities. For the scalar modes, we obtain conditions 
for avoiding ghost degrees. In the first branch, we find that one of the scalar modes becomes 
a ghost at the late stages of the evolution. Conversely, this problem can be avoided in the 
second branch. However, we also find that the constraint for the second branch prevents 
the doubly coupled matter fields from being the standard ingredients of cosmology. We thus 
conclude that a realistic and stable cosmological model requires additional minimally coupled 
matter helds. 
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1 Introduction 

Over the last century cosmology has progressively become an empirical scientihc discipline. 
With the high precision achieved by the cosmological observations, cosmology is now consid¬ 
ered as an indispensable arena to test fundamental physics. Nobel prizes have been awarded 
to the astonishing discoveries by the measurements of supernovae and Cosmic Microwave 
Background (CMB) radiation, which consolidated the standard model of cosmology. Another 
irreplaceable cosmological information is provided by galaxy surveys which is used in CMB 
secondary anisotropies and weak gravitational lensing. ESA’s PLANCK satellite releases very 
high precision measurements of the CMB and EUCLID mission will provide overwhelming 
amounts of data with very good control of systematical errors [1-3]. All these probes allow us 
to disentangle the composition, structure and evolution of the Universe. Even if the standard 
model of cosmology accounts for most of the cosmological observations, there exist anomalies 
in tension to this model [4] and theoretical discontent about the unnatural smallness of the 
observed cosmological constant [5]. 

The aforementioned difficulties have initiated a major motivation to study modihcations 
of gravity in the IR with additional dynamical helds. One of the most natural ways to modify 
gravity in the IR is to give the graviton a mass. In this context the de Rham-Gabadadze- 
Tolley (dRGT) theory has received much attention, being the hrst formulation of a non-linear 
covariant theory of massive gravity with the correct degrees of freedom [6, 7]. The potential is 
constructed in a way that guarantees the absence of the extra degree of freedom, dubbed the 
Boulware-Deser (BD) ghost [8]. The theory can be promoted to a massive bimetric theory 
by the inclusion of the kinetic term for the reference metric [9]. After constructing consistent 
self-interacting spin-2 theories as in bigravity, a mandatory question to address is how these 
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fields couple to the matter fields [10-24], In bi-gravity, both metrics are put on an equal foot¬ 
ing, and coupling matter to both metrics simultaneously might thus appear natural at hrst 
sight. However, as shown in [15] minimally coupling the matter sector to more than one metric 
generically introduces the propagation of ghost-like degrees of freedom already at the classical 
level. Moreover, quantum corrections at one loop destroy the special structure of the potential 
at an unacceptable low scale. A consistent way of coupling the matter sector to the two metric 
is to minimally couple the matter helds to just one metric, but not to both simultaneously. 
If one insists on coupling the matter sector to both metrics, a new composite effective metric 
constructed out of both metric was proposed in [15]. A nice feature of the coupling is that 
quantum corrections at one loop maintain the nice potential structure. Furthermore, this new 
matter coupling evades the no-go result for the flat Friedmann-Lemaitre-Robertson-Walker 
(FLRW) background present in the original formulation of massive gravity and also allows the 
propagation of the hve physical degrees of freedom of the graviton sector without ghost and 
gradient instabilities [25]. From the lessons learned about the quantum corrections coming 
from matter loops, one can construct other classes of new effective metrics through which the 
matter held could couple to both metrics at the same time [23]. 

The absence of the BD ghost strongly relies on the relative tuning of the potential interac¬ 
tions. It is an unavoidable question to address whether or not this relative tuning is radiatively 
stable. These questions were addressed in [26, 27] for massive gravity. In the decoupling limit 
the theory is protected from quantum corrections thanks to the non-renormalization theorem 

[26] (in the same way as the Galileon interactions are protected from quantum corrections 
due to a similar non-renormalization theorem [26, 28-31]). The antisymmetric structure of 
the potential interactions in the decoupling limit guarantees that any vertex will contribute 
at least with two additional momenta applied on the external leg to the transition amplitude 
and therefore the classical interactions remain untouched. Beyond the decoupling limit, the 
graviton loops do detune the specihc structure of the potential interactions, however the mass 
of the introduced BD ghost is never below the cut-off scale of the theory and hence is harmless 

[27] . The same detuning happens in the bimetric generalization of the theory [23]. 

The phenomenological aspects of the dRGT massive gravity have been widely studied, 
specially the potential impact on cosmology [32-52], In the decoupling limit of the theory, 
it has been shown that the theory admits self-accelerating solutions [32], What essentially 
happens is that the helicity-0 degree of freedom of the massive graviton forms a condensate 
whose energy density sources self-acceleration. Furthermore, small fluctuations around these 
self-accelerating backgrounds are stable. Another crucial result is that the fluctuations of the 
helicity-O held do not couple to the huctuations of the helicity-2 held (hence to the matter 
helds since it is the huctuations of the hehcity-2 held which couple to the matter held). Thus 
the cosmological evolution is exactly as in the standard AGDM model without the need of 
the Vainshtein mechanism [32], Nevertheless, these self-accelerating solutions found in the 
decoupling limit do not guarantee the existence of full solutions with identical properties in 
the full theory. One can consider the solutions in the decoupling limit just as a transient state 
of the full solution. 
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For the full theory, the construction of a realistic and stable cosmology has been a chal¬ 
lenge so far. For dRGT massive gravity with Minkowski reference metric, there are no flat 
FLRW solution [33] although self-accelerating open FLRW solutions exist [53]. However, these 
self-accelerating solutions (and similar ones with other reference metrics) are known to have 
three instantaneous modes [39] and lead to a nonlinear ghost instability [34]. Although for a 
de Sitter reference metric, other solutions have also been found [43, 54], these suffer from a 
Higuchi [55] type ghost at high energies [43]. The attempts to circumvent these issues roughly 
fall into two categories: i) breaking the FLRW symmetries in the hducial metric which do not 
affect the FLRW symmetries of the physical metric [36, 37, 40, 42, 56-58], or breaking these 
symmetries completely while relying on a cosmological screening mechanism [33] ; ii) extending 
the theory by adding new degrees of freedom [59-61]. 

The bimetric theory with dRGT tuning proposed by Ref. [9] can be seen to fall into the 
second category above. The cosmology of the bigravity theory has been already studied in 
a multitude of very interesting works [11, 12, 20, 58, 62-72]. However, just like in massive 
gravity, attaining a stable cosmology is highly challenging. For minimally coupled matter 
helds and a small interaction term m ~ Hq, the theory admits several branches of solutions: 
a self-accelerating branch which has three instantenous modes and is thus unstable [73], a 
branch where there is an early time gradient instability [73] and a branch where there is a 
crossing of a curvature singularity [74] (dubbed “inhnite” branch in Refs. [70, 72]). In the latter 
branch, even if one neglects the consequences of the singularity, three degrees (two tensor and 
one scalar) become ghosts at the early stages of the evolution [75, 76]. On the other hand, 
for a strongly interacting bimetric theory m 2> Hq, a viable but hnely tuned solution exists 
[75, 77]. 

In this work we will study the cosmological implications of the new effective metric pro¬ 
posed in [15] in the context of massive bigravity. The background evolution for one branch of 
solutions was already studied in [20] and we will push further this analysis to the level of the 
perturbations, while also paying attention to the other existing branches of solutions. After 
quickly reviewing the bimetric generalization of dRGT theory and clarifying our conventions 
and notations in Section 2, we hrst study the background equations of motion on FLRW 
space-time in Section 3. We comment on the existence of two branches of solutions and their 
implications focusing on their early and late-time behaviors. In the same way as in massive 
gravity [25], the new coupling enables us to avoid the no-go theorem for FLRW solutions. In 
Section 4 we turn our attention to the stability of the perturbations and show that even if 
the tensor perturbations can be maintained free of ghost and gradient instabilities, the vector 
and scalar perturbations unavoidably yield gradient and ghost instabilities at early and late 
time evolutions respectively in one of the branches of solutions. We then conclude in Section 5. 

Throughout the paper, we will denote traces by [...], for example the contractions of 
rank-2 tensors as = [/C], = [/C^] = (/C^i,)^, = [JC^] = (/C^,.)^ etc. 
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2 Massive bimetric gravity 


We will start this section with a review of the massive bimetric theory and collect all the 
important quantities. Our action is given by 

^5 = / -^y/^(R[g] + 2m‘^^anU[}C]j+y/^^Cx{ges,x) 


2 

j -^\/~/-R[/] + \/—f RraaXteAf] + '^matter[i?] 




( 2 . 1 ) 


with the very specihc potential interactions [6, 7] 

Mo|K) = = 1 

Ui\K] = = [K] 

U^K] = = 1 (1K]= - [K?]) , 

«3|K| = = 1 (1K]= - 3|K|1K"] + 2|K=>]), 

i^4|K| = i (1^1“ - 61K]=1K=] + SIK^ + 8lK]|K=>| - 6lK‘]) 

( 2 , 2 ) 


where the square brackets denote trace operation and 8 is the Levi-Civita tensor. The matrix 
/C is given by 

Keij,/] = «-(v<rv )7 (2.3) 

In [15] a new effective composite coupling was proposed 

gfu = + 2aPgaf, ^ ■ (2.4) 

In this paper we will consider a scalar held y with a generic kinetic term of the form 

C^ = V^PiX^), (2.5) 


where stands for 


= -9es9f,xduX ■ (2.6) 

The energy density, the pressure and the sound speed of the scalar held correspond to 

P(Xy)' 


Px = 2/’(AX)'.'f, - P(XJ, P, = P(XJ, 


7 2F(Jf,)"AX + F(A,)' ’ 


(2,7) 
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where prime denotes derivative with respect to the argument. Instead of a k-essence held [78] 
as in (2.5), one could also consider a more general scalar held, for instance a Galileon held [79] 
or a Horndeski held [80] or even a general vector held [81-83], but just for simplicity we will 
restrict ourselves to a k-essence held. In our setup we will consider two cosmological constants 
as a place-holder for the matter helds which either couple only to g or to 

■^matter = ~^g^g i -^matter = ~^f^f ' (2-8) 

In the following section we will hrst study the background equations and discuss the diherent 
branches of solutions together with their early and late time regime behaviors before moving 
on to the perturbations. 

3 Background evolution 

The bimetric theory in the original formulation without the composite coupling has been 
already intensively studied in the literature. In the context of the new ehective coupling, 
Enander, Solomon, Akrami and Mortsell have found cosmological solutions that resembles 
ACDM for a given choice of parameters [20]. This opens up new and exciting avenues to 
explore. In this work, we will perform perturbations about these cosmological solutions and 
discuss their stability. We assume homogeneous and isotropic flat FLRW metrics 

ds"^ = -|- a^gdijdx^dx^ , 

dsj = —n^fdt^ -\- Ojdijdx^dx^ . (3.1) 

Whit our Ansatz the effective metric becomes 

dsls = -Nl^dt^ + alg6ijdx’‘dx^ , (3.2) 

where the effective lapse Ngg and scale factor Oefr are dehned as 

Nes = aUg + (3nf , a^s = aag + (3 aj . (3.3) 

Compatible with our above homogeneous and isotropic Ansatz the y held also only depends 
on time x = xit)- The introduction of the three Hubble rates 

^9 — ^g/i'^g^g)^ = 'tz (8-4) 

JVeff Oeff 

and the ratios of the background scale factors and lapse functions 

r = (3.5) 

ttg '^gO'f 

^ At this stage, we would like to point out that the introduced cosmological constants Ag and A f serve only 
as a technical tool to keep track of the energy densities of additional fields and in no way should be thought 
as a real representation of physical matter fields. It would be important to study cosmological solutions and 
perturbations on them taking more realistic matter fields into account. 
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and the following function encoding the potential interactions 

U{A) = —ao + A{A — l)ai — 6(^4 — 1)^02 + 4(^4 — 1)^Q!3 — {A — 1)‘^«4 (3.6) 

will be benehcial for the exposure of the equations. Throughout this paper, we assume A > 0 
and r > 0. The action (2.1) in the mini-superspace becomes: 


S 

— = M„ 
V ^ 


^ ' dta^Ug 


dt alffNesP 


3Hg — m^p, 


m,g 


Mi 


dt 


a/Hf 


Af-3Hi 


m 


2^ 

Mj 




PmJ 


( 3 . 7 ) 


where for convenience we used the following functions as effective energy densities from the 
mass term instead of the parameters an 


Pm,giA) — U{A) ^ 1 Pmj{A) — ^^3 


( 3 , 8 ) 


We can now compute the background equations by varying the action (3.7) with respect to 
Uf, Qg, Qf and y. The variation with respect to the two lapses yield the Friedmann equations 


3 Hg — Ag Pm,g + 


a a 


eff 


M>3 


Px 


3 Hj = Af + + 


/3 a; 


eff 


a3 713 


Px ' 


(3,9) 


On the other hand, varying the mini-superspace action (3.7) with respect to the scale factors 
ttg and Qf, then using the Friedmann equations results in the following acceleration equations 


2 Hn 


rin 


2H^ 


= rn? J A{r — \) 


a a 


eff 


M>3 ^ 


, Nes/ Oeff n 

Px i - 1 - Px 

ng/ag 


m^M^J Ida 

(r-1)- 


3 

eff 


713 r 


Mja^A^ 


1 Nes /Oeff 

Px i -^- Px 

r Ug/ttg 


(3.10) 


where we have dehned the function J = \dApm,g for clarity. Last but not least, the equation of 
motion for the y held is just the standard conservation equation for a held minimally coupled 
to the gfeff metric 


1 1 d 


X 


Wff dt ViVeff 


+ 3i7eff 


X 

N, 


eff 


{px + Px) 


(3.11) 


Combining the equations for the lapse Ug, the acceleration equation for ag and the equation 
of motion for the scalar held x yields the following constraint equation 

m^JAH,) = Q. ( 3 . 12 ) 




From this constraint we immediately see that there are two branches of solutions. 
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3.1 Branch I 

The first branch of solutions corresponds to the case in which 


H, -AH, = 0 


(3,13) 


in the constraint equation (3.12). Combining the equations of motion for the lapses in (3.9) 
gives rise to the algebraic relation 


m 


2 





+ A, - K,A^ 


Px 

Ml 




(3.14) 


Another additional constraint comes from ^ (Hg — AHj) = 0. Using the acceleration equa¬ 
tions (3.10), this constraint can be equally written as 


2(r - 1)W 


{aes/agf ( f l3M^ \ N^s/aes p 

M2 MjAj Hg/ag ^ 


(3.15) 


where we have defined 


W = 


m^(M^ + M^A^)J 

2M|A 



(3.16) 


These two consistency relations ensure that the metric evolves according to the constraint 
Hf = H/A, by fixing A and r. 

The relation (3.14) can be used to analyze the early and late time limits of the cosmo¬ 
logical evolution. Assuming no bare cosmological constant (A^ = Aj = 0) and that at late 
times, the contribution from the two-metric interaction term cx dominates the expansion, 
we require 

^—:——T Pyx <^1, late times . (3.17) 

al im?M^ 1 V / 

This requirement can be consistent, provided that the energy density of the matter held y 
decreases with expansion. For instance, for pressureless dust, Py. oc Oeff- If the late time 
acceleration purely stems from the m‘^Pm,g term, i.e. parametrically ~ Hq, this means that 
the contribution from py. is already negligible at the present time. In this late time regime, 
the consistency relation (3.14) then gives 


m 


2 








0 . 


(3.18) 


or, at late times, A approaches to a constant value determined by the root of the above 
equation. Then equation (3.15) implies that r ~ 1. This asymptotic regime is self-accelerating. 
Similarly, at early times, one has 

3 

Px » 1 > early times. (3.19) 
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If for the sake of argument, we assume no other matter held minimally coupled to a single 
metric, the consistency relation (3.14) implies that a — ~ 0 [20]. 

Of course, these constraints will be drastically modified if there is any “regular” matter 
sector, in the sense that it couples minimally to only one of the two metrics. In this case, 
the right-hand side of (3.14) may be dominated by these matter helds rather than the mixed 
coupled one, reducing to [75]. 


3.2 Branch II 

The second branch of solutions corresponds to the case where 


m^J — 


=0 

My, ^ 


(3.20) 


coming from the constraint equation (3.12). The dynamics of this branch is drastically different 
than the one for usual minimal matter couplings (a = 0 or /3 = 0), in the sense that the matter 
sector directly enters the constraint. On the other hand, for an arbitrary matter held coupling 
and non-zero pressure, the solutions A are generically time dependent. In this branch, the 
acceleration equations (3.10) yield an interesting consistency relation 


H f _ 

Ifg ~ MjaA^' 


(3.21) 


Also, requiring that the constraint (3.20) is satished throughout the evolution, i.e. — 

a(3{aes/{Mgag))‘^P^'j = 0 gives another constraint equation in this branch 


2171"^ a gUg 
(r — 1)A 


[-rAJ{Hg 


where 


AH,) + T{H,- tAH,)] + [2P^ - 34(p^ + FJ] = 0, 

(3.22) 

T = AJ+^^^M1^^Ja. (3.23) 


Let us briehy discuss the early and late time limits in this branch. Similarly to the 
arguments for the branch-I solutions we made above, we assume A, = Af = 0. Then, at late 
times we assume that and decays to values sufficiently smaller than nPMg:'^ 


nPM'^a^g 




nPMgO^g 


Px < 1 • 


(3.24) 


^ We remark that this assumption is not the unique one that should hold in the late time, and there could 
be another late time regime in which these conditions are not satisfied. In fact, when parameters at are such 
that J{A) > 0 for any T, there is a minimum value of Py, which can parametrically be of order Mgim?. In 
this case, there is no regime where the y is neither dominant nor sub-dominant. As we expand in the next 
Section, these backgrounds can evade the strong-coupling problem we will encounter. 










In this regime, equation (3.20) implies J 1, and then A converges to the constant values 
dehned J = \dApm,g = 0, that is, 


^ _ (^2 + 2^3 + «4 i •\/q' 2 + a3(Q!3 — Q!i + Q!2) ~ (o^l + tt2)ci4 

tt3 + tt4 

As a result, the effective energy densities pm,g and pmj are also constant, leading to a late 
time self-acceleration terms in the Friedmann equations for both metrics (3.9). Then, if we 
further assume that A, r = 0(1), the constraint (3.22) reduces to 

(H, - tAH,) ~ 140,4 (H, - vAH,) « 1 , ( 3 . 26 ) 

which suggests that Hg ~ rAHf in the late time regime in a generic case where = 0(1). 
This can be also easily seen by the relation 


A 

HgA 


H/rA - Hg 


(3.27) 


which for a constant A yields the said relation at late times. 

In dRGT massive gravity and bimetric theory, this behavior with constant A (3.25) 
is valid at any stage of the evolution since the constraint is unaffected by the minimally 
coupled matter fields. Although in this case, this branch is known to contain three degrees of 
freedom with vanishing kinetic terms [39, 69, 73], we will show in the following Section that 
for a composite matter coupling, this problem can be evaded even at the late stages of the 
evolution where the situation is similar to standard dRGT. 

The early time regime could be dehned by the condition that the x field is the dominant 
component of the evolution, i.e.. 


Q/^Qeff 


P. » 1, 


<^f3als 

nHMga^ 


Px > 1 • 


(3.28) 


To see if these conditions could be satished, it is useful to rewrite the constraint equation 
(3.20) into 


MgTfi cxi -\- 3 q (2 + 3 q (3 -|- 0(4 — 2( q(2 + 2 q(3 -|- a^jA + ( 0(3 + 0 ( 4 )^^ 

(« + /3A)2 • 

We may view this equation as an algebraic equation to determine A for a given P^. Since 
we assumed A > 0, the right-hand side is bounded and thus solutions to this equation exist 
only when parametrically |P^| < m?Mg if Q(,/?, Q(o,i,2,3,4 = C>(1) During the course of the 
evolution, the pressure of the y matter held reaches a maximum value. At this extremal value 
of the pressure, we hnd from Eq. (3.22) that 


A 

ria 


(a 


Am, 


P^=P^(max) 




(3.30) 


^ Here, we suppose that cx/? > 0, which is necessary for the stability of the vector and scalar perturbations 
as we shall see in (4.13) (with (3.20) imposed) and (4.20). 
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which further implies 


P^=P^(max) 


0. In other words, the effective metric bounces at 


this instant. As a result of this bounce, the x matter field never dominates the expansion. 
This result suggests that one cannot have a realistic cosmological model in this branch, unless 
“regular” matter fields coupling minimally to either of the two metrics exist. When such 
matter fields are included, the bound given by Eq. (3.29) will not change, but it will still be 
possible to have a matter or radiation dominated stages using the added fields. 


4 Perturbations 

The previous section was dedicated to the study of the background equations of motion and the 
existence of flat FLRW solutions in the presence of the new coupling to the matter sector in the 
framework of bigravity. This new coupling had tremendous implications in the framework of 
massive gravity since the original formulation does not accommodate for flat FLRW solutions. 
It was shown that the new coupling can circumvent this no-go result and that the BD ghost is 
not present around these flat FLRW solutions [15]. Additionally, the decoupling limit analysis 
had shown that the theory is free of the BD ghost at least up to the strong coupling scale 
Ag = Mgrri^ and hence provides a valid effective field theory up to this scale. Furthermore, the 
perturbations around these cosmological solutions were investigated in [25] and the stability 
conditions for all the propagating six degrees of freedom were explored. The theory also avoids 
the no-go result of vanishing kinetic terms for the physical degrees of freedom in the case of 
the standard coupling. In this work we will have a similar analysis as in [25] but for the case 
of bigravity. In order to determine the stability conditions around the flat FLRW solutions, 
we consider the following Ansatz for the two dynamical fields 

Qgudx^dx’^ = + 2$) dt^ + 2ngagVi dtdx^ + o?{5ij + hij)dx''dx^ , 

fg^ydx^dx'^ = — nj{l + 2(p) dt^ + 2nfafVi dt dx^ -|- -|- 'jij)dx^dx^ , (4.1) 

where the perturbations (<f), V), hij) and (</?, Uj, 7 jj) are functions of time and space. We will 
further decompose the tensor perturbations into their trace and traceless parts 

1 1 

hij — —dijh - 1 - ) ( 4 - 2 ) 

Furthermore, we consider the following perturbation for the matter sector y 

X = Xo{t) + Sx. (4.3) 

It will be convenient to decompose the perturbations in Fourier modes with respect to the 
spatial coordinates 

/ d^k 

j^;^Qj:{t)exp{ik-x) + c.c., (4.4) 

where Q represents the perturbations (<h, Rj, hij) and ((p, n*, Xij) respectively. We will perform 
the stability analysis of the perturbations for each sector separately and focus on the conditions 
arising from the absence of ghost and/or gradient instabilities. 
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4.1 Tensor perturbations 

First of all, we will start with the analysis of the transverse-traceless part of the metric 
fluctuations. Our Lagrangian (2.1) together with the Ansatz (4.1) for the metric perturbations 
in Fourier modes becomes 


5, 


( 2 ) _ 


9 


tensor 


d^kdiUga^g 


h* 

ij,k k 

rii 


e Mj r 

-h*. rh^d + d 


a 


2 ij,k k 


M 2 

9 


-mis 


■ ★ * ^7 

'y -*'y-dd 
ij,k k 






ij,k 'ij,k. 


, (4.5) 


where mis stands for the shortcut notation 


mis = 


= m^F — a[3A ( — 

VOg 


A^eff/ ®eff Px 

Ug/ttg 


M 2 

9 


= m^T 


a(3A{a + (3A){a + (3rA)^^ , (4.6) 

M^g 


and hij and 7 *^ are assumed to be transverse-traceless in this section. As one can see in the 
expression of the quadratic action (4.5), the tensor perturbations have the right sign for the 
kinetic term. The tensor modes do not impose any further constraint, since r > 0 and A > 0 
already. They also do not give rise to any gradient instabilities. In the late time regime, 

the quantity —j^^^jj^mls can be identihed as the effective graviton mass squared [75]. We 

discuss its properties in the next section. Even though the tensor perturbations are free of 
ghost and gradient instabilities, there can be tachyonic instabilities when mis < 0- However, 
the parameters of the model can be chosen such that the timescale of the instability would be 
harmless. For more details on the tachyonic instability please see Ref. [84]. 


4.2 Vector perturbations 

We now concentrate on the vector perturbations. We decompose the metric perturbations as 


h-■ t = - (kiE ■ r + kjE-^ , 

ij,k 2 \ 9 i,kJ 1 


where d^Ei = d^Bi = d^'bi = <9*S'* = 0. After taking into account the background equations of 
motion, the quadratic action of the vector perturbations becomes 


5, 


( 2 ) 


M2 


vector 


d^kdt Ugtt^g 



2BA 

i.k 


+ 





(A 

2 M 2 \ 


«/ y 

V 


^mBik- 

’■Kf) . 


2hA 

i,k 

«/ 


k'^m^ 


eff 




■k 

i,k 


^*k) 

i.k' 


(4.8) 
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where we have introduced the following dehnition for convenience 



r + 1 ttiIq 


m^J + ajS 



( 1 -^eff/Oeff Px Px I Px A 

Vr + 1 n,/a, M^J 


(4.9) 


This quantity is proportional to the squared sound speed as discussed below. Naively counted 
there are four vector modes (i?i, hi, Si and Ei). However, the vector modes Bi and hi are 
not dynamical degrees of freedom and thus we can use their equations of motion in order to 
integrate them out. By doing so we obtain for the non-dynamical degrees of freedom 




2ng 2ng 


where € is given by 


(Kk - 






r _ i,k . 


^ -(E. 

' V 1 . 


S. 






^9 Ke 


+ 


Ml 


MjA^ 


r + 1 


-1 


. (4.10) 


(4.11) 


Integrating out the non-dynamical vector modes Bi and hi through equations (4.10) results in 


c(2) 

^vector 



(fkdtngO?g 



1 + MjA^/Mj 


m, 


eff 



where we dehned the gauge invariant combination 





i.k 




(4.12) 


For the absence of ghost instabilities we have to impose that the kinetic term has the right 
sign. As one can see from the expression for the quadratic action of the vector perturbations, 
this is guaranteed by the requirement C > 0. A sufficient condition for this to hold is 



r -|- 1 
2A 


m^J + a(3 



f 1 A^eff/o-eff Px A Px ^ 

Vr + 1 Ug/ag M2 M^J 


>0 , 


(4.13) 


which is satisfied if J > 0, > 0 and p^ + Py^ > 0. Also the squared sound speed for the vector 

mode is given by should be positive to avoid the gradient instability. For 

this stability, on top of the conditions to avoid the ghost instability, another condition rrigg > 0 
needs to be satished. 

Let us briefly study the late time limit discussed in section 3. In this limit, mlg behaves 
for both solution branches as 

mlg ~ m^F , (4.14) 

and Cy behaves as 



r 

AJ 


m+l+r)+ 

2aPA[l+r+P^^) 


J Px~^^x 


(Branch I) 
(Branch II) 


(4.16) 
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For the branch I, rUes and Cy converges to constants for generic choice of the graviton mass 
parameters Q!j=o,i, 2 , 3 , 4 , and also gradient instability is absent if F > 0 in addition to the no¬ 
ghost conditions is satished. The same statement applies also to the branch II, that is, the 
gradient instability is avoided if F > 0. 

For the branch II, the sound speed generically becomes divergent at late times as Cy oc 
{py- -|- Px)~^, indicating the presence of strong coupling and superluminality in this model. 
This is reminiscent of the vanishing kinetic terms in the minimally coupled case; the late time 
limit assumed in (3.24) effectively reintroduces the very problems we wish to avoid. Such a 
divergence can be avoided if F converges to zero at late time, which may be realized if the 
parameters at are hne-tuned. A drawback of such a hne-tuning is that the effective graviton 
mass mlg ~ m^F becomes zero and then the effect of the graviton mass would be negligible 
at late time. Alternatively, further tuning of the evolution such that r/A^ S> 1 would result 
in a large mass gap larger than the strong coupling scale in the vector sector, as we can 
see from Eq. (4.11). Finally, one can also choose the parameters such that p^, + never 
becomes too small. In other words, it is possible to discard the assumption that a regime 
(3.24) exists. For example, such a situation is realized if we choose the parameters such that 
-|- 3(^2 -l- 303 -|- 04 > 0, 02 -l- 2 a 3 -|- 04 < 0 and 03 -|- 04 > 0. As a result J{A) > 0 for 
any value of A and through the constraint (3.20), P^ is bounded from below parametrically 
at order rrigMg. As a result, the x matter never becomes subdominant over the two-metric 
interaction term. In the massive gravity version of this theory, the cosmological solutions 
studied in Ref. [25] have a similar behavior. Both py. and P^ are constrained to be of order 
nPMg and generically there is no y-held domination nor subdomination over the mass term. 

In the early time regime for the branch I, the no-ghost condition is kept satished as long 
as J > 0 is maintained. However, and then Cy become negative when is sufficiently 
large, as we can see in equation (4.6) in which F stays hnite since A converges to a constant 
in the early time regime. Thus there would be a gradient instability in the early time regime 
for the branch I. 

4.3 Scalar perturbations 

As next, we will concentrate on the stability of the scalar perturbations. The metric pertur¬ 
bations can be decomposed into scalar perturbations as 

~ 'i'kiBj :, g — i^ikj —k ^ , 

^^,k = ikih > 7p-fc = - (^kikj - ^k^^ . (4.16) 

As one can see nine degrees of freedom appear in form of a scalar held <F, cp, B, b, E, S, S, ^|J 
and 6x from which only two are dynamical. In order to see this we will integrate out the 
non-dynamical degrees of freedom and also hx the gauge. First of all, let us introduce the 
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following gauge invariant variables [75] 
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(4.17) 


We have the freedom to choose a gauge and remove two degrees of freedom. We choose the 
gauge hxing in a a-/? symmetric way as follows 

I3A 


Ej^ = 


P 

a + 13 A ^ ’ 


Sk = -^A^ ^k = 


ol (3 A 


ol f3A 


V'fc, Sr = 


a 


(X f3A 


% ■ (4.18) 


Under this gauge condition, the dynamical variables reduce to (Yi,y 2 ,h 3 ) = {Sx,E,'ijj). We 
can further integrate out the non-dynamical degrees of freedom $, 0, i? and b. After using their 
equations of motion, we hnd that is missing the kinetic term and is non-dynamical. It is 
identified as the degree of freedom of the would-be BD ghost. The cumbersome mathematical 
expressions involved make it difficult to present it in a manageable way. In the following, we 
will concentrate on the UV limit and for simplicity set = Aj = 0 without loss of generality. 
The kinetic term in the UV limit after integrating out the would-be BD mode becomes 


c(2) 

^scalar 


9 



dtd/kug ttg 


Y- Y3' 

-AJC^ 


n„ 


n„ 


(4.19) 


where the components of the kinetic matrix are given by 


/Cii = 


1 als Ug 

clX^ al N,s 


M2 
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/Cl2 = 0{k^) 

^ ^ + ft) - + 1)] 


-1 


+ 0{k^ 


(4.20) 


In the UV regime, we find that the kinetic matrix is diagonal at leading order. Therefore, the 
conditions for avoiding the ghost instability are given by /Cn > 0 and /C 22 > 0. The latter one 
is equivalent to 

Aa^Acl^alsiP^ + p^) - M^almls{r + 1) > 0 , (4.21) 

if the stability condition for the vector perturbations Cy > 0 is imposed. 

Let us examine the no-ghost condition in each solution branch. In both branches, the 
condition /Cn > 0 is satished if the energy condition -|- P,,. > 0 is satished. As for the other 
condition /C 22 > 0, in the late time regime for the branch I it reduces to 


/C22^ 


k\P^ + p^y 


2M4 


(r -I- I)2a2m2 J 


< 0 


(4.22) 
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which indicates that this branch is plagued by the ghost instability in the late time regime. 
On the other hand, in the early time regime where Py. and are much larger than nPMg, 
the kinetic term /C 22 reduces to 


/C22^ 


k^a(5Aal^{p^ + Pxf 
AM^ 

k^al3Aalg{p^ + 

AM^ 


1 + r 


^ sff/ ^efF ^ 


Ug/ag J 

r + MjA^/Mf 


{Px + ^x) + (1 + ^)Px 


-1 


1 + M^‘^/M^ 


(Px + ^x) + (1 + ^)-Px 


-1 


> 0 


(4.23) 


where we used a ~ that holds in this regime to show the latter equality. Thus the 
branch-I solutions are free from the ghost instability in the early time regime. 

For the branch-II solutions, using equation (3.20) we can show that K .22 reduces to 


/C22 — 


k‘^{Py- + p^) a(]Aa\ 


eff 


4M2 


ag{ar + (]A) 


> 0 


(4.24) 


and thus the branch II is free from the ghost instability if the energy condition > 0 

is satished. Although we have only obtained the large momentum expressions for the kinetic 
terms, one can expect that in principle, this branch can lead to a strong coupling at late times. 
However, as discussed in the previous Subsection, one can either generate a mass gap larger 
than the strong coupling scale, or prevent the x field to ever become subdominant against the 
two-metric interaction term, by carefully choosing the parameters of the theory. 


5 Conclusions 

We have studied the cosmological perturbations in the massive bimetric theory with the 
new matter coupling proposed in [15]. The background equations give rise to a constraint 
(3.12) which divides the solutions into two branches. We summarized the late time behaviors 
of the two solution branches, while the early time behaviors were given with the cautious 
assumption that the doubly coupled matter y is the only matter held in the model. The 
branch-I solutions dehned by Hg — AHf = 0 has been studied by [20]. The second branch 

dehned by — j§^Px = 0 is typically dismissed on the basis that various problems 
including strong coupling and nonlinear instability may be present. However, we hnd that 
this branch may give a well-behaved cosmology at late times; at early times however, as 
the doubly coupled matter is not allowed to dominate, the standard model helds need to be 
coupled minimally to one of the metrics. 

In the latter part of this work, we conducted quadratic cosmological perturbations and 
clarihed the conditions to avoid the ghost and gradient instabilities for a general background, 
and then analyzed those conditions for each solution branch focusing on the early and late 
time regime. We hnd that the tensor perturbations are free from both ghost and gradient 
instabilities if the energy conditions {p^ > 0, p^ + P^ > 0) and some additional conditions 
(T > 0, J > 0) are satished, while the vector sector for the branch I may suher from the 
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gradient instability in the early time regime. Abont the scalar sector, we analyzed the no¬ 
ghost conditions and fonnd that the energy conditions gnarantee the positivity of one of the 
kinetic terms, while the second scalar mode becomes a ghost in the late time regime of the 
branch-1 solntions. 

In the present work we have assnmed a single matter held which conpled to both metric 
simnltaneonsly. In the presence of additional matter minimally conpled to a single metric, 
the backgronnd solntions at early times wonld be modihed. Generalization to snch a case 
with additional matter helds wonld be an interesting target for farther stndies. The late 
time regime in the present stndy however provides a good approximation for some of the 
degrees of freedom of snch a complete model, since the backgronnd evolntion is dominantly 
governed by the mass term and the effect of the additional matter helds wonld be small. 
Also, dne to complexity of the non-diagonal and time dependent action, we were not able to 
fnlly diagonalize the scalar sector. Therefore, we cannot present the freqnency eigenstates nor 
address the qnestion of gradient stability for these modes. It wonld be important to clarify 
these conditions to examine viability of the branch-II solntion which is free from the ghost 
instability. 

Note added: While this work was being completed, we became aware of Ref. [84], which 
also stndied the pertnrbative stability of the bimetric theory with donbly conpled matter. The 
pertnrbation analysis in Ref. [84] is performed at the level of the eqnation of motion, while 
the present stndy has been performed directly at the level of the action. The results are in 
agreement, while ghost instability in the late time regime of the branch-I solutions is shown 
for the hrst time in this work. 
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